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ABSTRACT 

We present ALMA observations of two moderate luminosity quasars at redshift 6. These quasars 
from the Canada-France High-z Quasar Survey (CFHQS) have black hole masses of ~ 10 8 M 0 . Both 
quasars are detected in the [C ii] line and dust continuum. Combining these data with our previous 
study of two similar CFHQS quasars we investigate the population properties. We show that z > 6 
quasars have a significantly lower far-infrared luminosity than bolometric-luminosity-matched samples 
at lower redshift, inferring a lower star formation rate, possibly correlated with the lower black hole 
masses at z = 6. The ratios of [Cn] to far-infrared luminosities in the CFHQS quasars are comparable 
with those of starbursts of similar star formation rate in the local universe. We determine values of 
velocity dispersion and dynamical mass for the quasar host galaxies based on the [C n] data. We find 
that there is no significant offset from the relations defined by nearby galaxies with similar black hole 
masses. There is however a marked increase in the scatter at z = 6, beyond the large observational 
uncertainties. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — quasars: 
general 


1. INTRODUCTION 

Improved astronomical observational facilities have en¬ 
abled the discovery and study of many galaxies at an 
early phase of the Universe’s history. It is now possi¬ 
ble to witness the majority of the stellar and black hole 
mass growth over cosmic time and identify how physical 
conditions at early times differ from now. One of the 
major relations to be determined as a function of time is 
the tight correlation between black hole mass and galaxy 
properties observed for nearby galaxies (see Kormendy 
& Ho| 2013 for a review). Observations of tms relation 
at high-redshift are critical to understanding the cause 
because most of the growth occurred at early times. 

Attempts to measure black hole and galaxy masses at 
high-redshift face a number of problems. Black hole mass 
measurements cannot be made directly by resolved kine¬ 
matics of gas or stars within the black hole’s sphere of 
influence, nor by reverberation mapping. Instead black 
hole masses, Mbh, of quasars can be measured at any 
redshift using the single-epoch virial mass estimator that 
involves measuring a low-ionization broad emission line, 
such as Mg II or H/3, and calibrating the location of the 
emi tting gas w ith low-z reverberation-mapped quasars 
(|Wandel et al.||1999|). For AGN with obscured broad 
lines Mbh can only be estimated from the luminosity 
making an assumption about the accretion rate relative 
to the Eddington limit. 

Measuring galaxy properties, such as luminosity or 
velocity dispersion, <r, of distant quasars is hampered 
by surface brightness dimming, the bright glare of the 
quasar and AGN (active galactic nuclei) emission line- 
contamination of spectral features. Up to z ss 1 there 
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has been considerable success in measuring AGN host 
galaxy luminositie s, morphologies and in some cases ve - 


locity dispersions (Cisternas et al. 
At higher redshifts (1 < z < 4) t. 


2011 Park et al. 2014). 
re galaxy light is more 


difficult to separate from the quasar, which, combined 
with greate r mass-to-light corrections, lead to larg er un¬ 
certainties ( Merloni et al.|20T0 Targett et al.|2012 ). The 
results of these studies are mixed with some evidence in 
favour of higher Mbh at a given galaxy mass. 

At yet higher redshifts it has pro ved impossible to mea- 
sure the galaxy light of quasars (Mechtley et al. 2012) 
before launch of the James Webb Space Telescope and 
instead the main method of determining galaxy mass is 


kinematics of cool gas in star-forming regions (Carilli & 
Walter 2013). Facilities such as the IRAM Plateau 3e 
Bure Interferometer, the Jansky Very Large Array and 
the Atacama Large Millimeter Array (ALMA) have suf¬ 
ficient sensitivity and resolution to resolve the gas in dis- 
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In partic¬ 
ular) ALMA has the sensitivity to probe z = 6 quasar 
hosts with star formation rates, SFR, in the tens of solar 
masses per year, rather than only in the extrem e star- 
bursts previously observable (Willott et al. 2013). The 
studies above focussed on z « 6 Sloan Digital Sky Survey 
(SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS) 
quasars with high UV and far-IR luminosities and found 
that their black holes are on average 10 times greater 
than the corresponding a for local galaxies, roughly con¬ 
sistent with a continuation of the evolution seen in lower 
redshift studies. 

Although observationally there appears to be an in¬ 
crease in Mbh with redshift at a given galaxy mass or 
(j, it has long been understood that there are selection 








































2 


Willott et al. 


biases that affect how closely the observations trace the 
underlying distribution. In particular, the steepness of 
the galaxy and dark matter mass functions combined 
with large scatter in their correlations with black hole 
mass mean that a high black-hole-mass-selected sample 
of quasars will have a systematic offs et in a towards lower 
values. This eff ect, fi rst identified by Willott et al. |2005|) 


and Fine et al. (2006) was studied in detail in Lauer et al. 
(2007) and numerous studies thereafter. The magnitude 
of tile effect depends upon the scatter in the correla¬ 
tion, which has not been conclusively measured at high- 


redshift, but app e ars to increase with redshift ( Schulze 
fe Wisotzkj||2014 ). Willott et al. (2005) and Lauer et al. 


(2007) showed that the bias is particularly strong for 
Mbh > 1O 9 M 0 quasars such as those in the SDSS at 
z ss 6 and therefore that the factor of 10 increase in Mbh 


(Walter et al. 

2004) could be accounted for by the bias 

(see also Schu 

ze & Wisotzki 2014). In comparison, there 


black hole masses of Mbh ~ 10 8 Mq ( Lauer et al.| [2007). 

An alternative to measuring the evolution of the as- 
sembled galaxy and black hole masses is to determine 
the rate at which mass growth is occurring. For quasars 
the bolometric luminosity is a measure of the black hole 
mass growth rate. For galaxies, the star formation rate is 
proportional to the stellar mass growth. The star forma¬ 
tion rate can be determined by the rest-frame far-infrared 
dust continuum luminosity. Additionally, the interstel¬ 
lar [Cn] far-infr ared emission line is we l l-correlated with 
star-f ormation (De Looze et al. 2014| Sargsyan et al. 


2014 1 so can also b e used as a star formation proxy. 

In" Willott et al. (2013, hereafter Wil3) we presented 
Cycle 0 ALMA observations in the [Cn] line and 1.2 mm 
continuum for two z = 6.4 quasars from the Canada- 


France-High-z Quasar Survey (CFHQS, Willott et al. 
2010b). These quasars have M B h ~ 10 8 M Q , a factor 
of 10-30 lower than most SDSS quasars known at these 
redshifts. One quasar was detected in line and continuum 
and the other remained undetected in these sensitive ob¬ 
servations placing an upper limit on its star formation 
rate of SFR< 4OM 0 yr -1 . 

In this paper we present ALMA observations of two 
further CFHQS quasars with similar redshift and black 
hole mass with the aim of providing a sample large 
enough to address the issue of how host galaxy prop¬ 
erties such as SFR, cr and dynamical mass depend upon 
black hole accretion rate and mass at a time just 1 billion 
years after the Big Bang. In particular, these quasars are 
not subject to the bias in the Mbh — cr relation discussed 
previously because of their moderate black hole masses. 
Cosmological parameters of Hp = 6 7.8 km s -1 Mpc -1 , 
Dm = 0.307 and = 0.693 (Planck Collaboration 2014) 
are assumed throughout. 


2. OBSERVATIONS 

CFHQS J005502+014618 (hereafter J0055+0146) and 
CFHQS J222901+145709 (hereafter J2229+1457) were 
observed with ALMA on the 28, 29 and 30 November 
2013 for Cycle 1 project 2012.1.00676.S. Between 22 and 
26 12 m diameter antennae were used. The typical long 
baselines were ~ 400 m providing similar spatial resolu¬ 
tion to our Cycle 0 observations. Observations of the sci¬ 
ence targets were interleaved with nearby phase calibra¬ 


tors, J0108+0135 and J2232+1143. The amplitude cali¬ 
brator was Neptune and the bandpass calibrators J2258- 
2758 and J2148+0657. Total on-source integration times 
were 4610 s for J0055+0146 and 5490 s for J2229+1457. 


The band 6 (1.3 mm) receivers were set to cover 
the frequency range of the redshifted [Cn] transition 
( I/ rest = 1900.5369 GHz) and sample the dust continuum. 
There are four ss 2 GHz basebands, two pairs of adjacent 
bands with all GHz gap in between. The channel width 
is 15.625MHz (17kms -1 ). 

The data were initially processed by North American 
ALMA Regional Center staff with the CASA software 
package 1 . On inspection of these data it became clear 
that the [C Ii] line of J0055+0146 was located right at the 
edge of the baseband, 1000 kms -1 from the targeted fre¬ 
quency defined by the broad, low-ion ization Mgn emis¬ 
sion line (zMgii = 5.983; Willott et al. 2010a). The Mgn 
line redshift is usually close to the systemic redshift as 
measured by narro w optical line s with a dispersion of 
270 kms -1 ( Richards et al.|2002 |. A large offset for this 
quasar was not particularly surprising for two reasons: 
firstly the signal-to-noise (SNR) of the Mgil detection 
is not very high and the line appears double-peaked due 
to noise and/or associated absorption; secondly the Lya 
redshift (z^ya = 6.02) is offset from Mgn by 1600kms” 1 
(in the same direction as the [C n] offset) and this would 
make the size of the Lya ionized near-zone negative, 
which is not physically sensible for a quasar with such 
a high ionizing flux and has not been obser ved in a sam¬ 
ple of 27 z ss 6 quasars (Carilli "et al.p OlO). Due to this 
redshift uncertainty the receiver basebands were set up 
so that the adjacent band covered the Lya redshift with 
zero gap between the two bands. 

The default ALMA Regional Center reduction ex¬ 
cluded 11 channels at each end of the 128 channel band. 
However, only the first 4 channels need to be excluded, 
so we re-reduced the data with CASA to include more 
spectral channels at the baseband edges. We checked 
that the noise does not increase in these extra channels, 
except for the very first and last channels to contain data 
so we excluded those. In summary our reduced product 
contains 118 of the original 128 channels per baseband, 
compared with 106 channels in the default reduction. 


3. RESULTS 

Figure |T] shows the reduced spectrum of J0055+0146 
from the two adjacent basebands. The final gap between 
the bands is only « 150 kms -1 and crucially the peak 
of the [C ii] line is contained within the higher frequency 
band. The lower frequency band contains only a small 
amount of the line flux but provides an important con¬ 
straint on the wings and hence the peak and width for 
a symmetric line. A single Gaussian plus flat continuum 
model was fit to the available data using a Markov-Chain 
Monte Carlo (MCMC). This process shows a good fit for 
a single Gaussian with FWHM = 359 ± 27 kms -1 . The 
formal uncertainty in the FWHM is very small consid¬ 
ering that there is some missing data. This is because 
a symmetric line model is used and with the peak and 
wings covered by data there is little margin for deviation 
in the missing channels. We add in quadrature an extra 
10% uncertainty in both the line flux and FWHM due to 

1 http://casa.nrao.edu 
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Fig. 1.— ALMA spectrum of J0055+0146 covering two adjacent 
basebands. The gap between the bands with no data is shaded 
in gray. The higher frequency band is centred on the redshift de¬ 
termined from the Mg II emission line whereas the lower frequency 
band covers the Lya redshift. The [Cn] line is found at the edge 
of the higher frequency band. The blue curve is a Gaussian plus 
continuum fit as described in the text. The red circle marks the 
continuum level independently measured in the three line-free base¬ 
bands. The upper axis is the velocity offset from the best-fit [C n] 
Gaussian peak. 


Velocity offset (km / s) 

1000 500 0 -500 



Observed Frequency (GHz) 

Fig. 2.— ALMA spectrum of J2229+1457 covering the single 
baseband containing the [Cn] line. The blue curve is a Gaussian 
plus continuum fit and the red circle the independent continuum 
level. The upper axis is the velocity offset from the best-fit [Cn] 
Gaussian peak. 

the missing channels. 

The spectrum of J2229+1457 is plotted in Figure [2} 
For this quasar the [Cii] line is centred in the band 
with no significant offset from the Mgn redshift. The 
line is consistent with a single Gaussian with a best 
fit FWHM = 351 ± 39 kms -1 , similar to the value for 
J0055+0146. The continuum level of the fit is again con¬ 
sistent with the independent continuum level determined 
from the three line-free basebands. Measurements from 
the spectra are given in Table [I] 

Figure [3] shows maps of the 1.2 mm continuum (white 
contours) plus the [C II] line (color scale) for the quasars. 
For both quasars there is no significant offset between 
any of the continuum or line centroids or the optical 
quasar position. The < 1" mm-optical offset is within 
the relative uncertainty of the optical astrometry. The 
more accurate [C 11 ] positions for the two quasars are 
00:55:02.92 +01.46.17.80 and 22:29:01.66 +14.57.08.30. 
For J2229+1457 there is only a marginal 2er detection 
of the continuum located coincident with the peak of 
the line emission. As seen in Figure [3] there are several 
other continuum peaks of this magnitude or greater in 



RA (J2000) RA (J2000) 


Fig. 3.— The color scale shows the integrated [Cii] line maps 
for the two quasars. White contours of 1.2 mm continuum emis¬ 
sion from the three line-free basebands are over-plotted at levels 
of 2,4, 6 cr beam -1 . The quasar optical positions are shown with a 
black plus symbol. The positional offsets between the optical and 
millimeter are most likely due to astrometric mismatch, rather than 
a physical offset. J0055+0146 is well-detected in both continuum 
and line emission. J2229+1457 has only a 2 a continuum detection 
that is spatially co-incident with the line emission. The restoring 
beam is shown in yellow in the lower-left corner. 


TABLE 1 


Millimeter data for the two CFHQS quasars 



CFHQS J0055+0146 

CFHQS J2229±1457 

^Mgll a 

5.983 ± 0.004 

6.152 ±0.003 

Z [CII] 

6.0060 ± 0.0008 

6.1517 ± 0.0005 

FWHM [CII] 

359 ±45 kms— 1 

351 d= 39 kms -1 

l[Cii] (Jykms -1 ) 

0.839 ±0.132 

0.582 ± 0.075 

^[CII] (Lq) 

(8.27 ± 1.30) x 10 s 

(5.96 ±0.77) x 10 s 

/l.2mm (jUJy) 

211 ± 34 

54 ±29 

£fir {Lq) 

(4.85 ±0.78) x 10 11 

(1.24 ± 0.67) x 10 11 

SFRjotq (Me yr- 1 ) 

83 ± 13 

60 ±8 

SFRfir ( M @ yr" 1 ) 

73 ± 12 

19 ± 10 

Gciii/^fir 

(1.70 ±0.38) x 10~ 3 

(4.80 ± 2.67) x 10" 3 


Notes.— 

a Derived from Mgn A2799 observations (|Willott et al.|2010a|. 
Uncertainties in Lfir, SFR[ CII ] and SFRfir only include measure- 
ment uncertainties, not the uncertainties in extrapolating from a 
monochromatic to integrated luminosity or that of the luminosity- 
SFR calibrations. 


the vicinity, so it is not considered a secure detection, 
however the measured flux and uncertainty are included 
in Table [I] 

The quasar 1.2 mm continuum flux-densities were con¬ 
verted to far-infrared luminosity, Lfir, assuming a typ- 
icalSED for high-redshift star-forming galaxies. As in 


Wil3 we adopt a greybody spectrum with dust temper¬ 


ature, T(j = 47 K and emissivity index, f3 = 1.6. To 
convert from far-IR luminosity to star formation rate we 
use the relation SFR (M 0 yr -1 ) = 1.5 x 1Q~ 10 Lfir fa ) 
appropriate for a Chabrier IMF (Carilli & Walter|2013|). 
We note that this assumes that all the dust contribut- 
ing to the 1.2 mm continuum is heated by hot stars and 
not by the quasar. An alternative estimate of the star 
formation rate comes from the [C III luminos ity. We 
adopt the relation in Sargsyan et al. (2014) of SFR 
(M 0 yr -1 ) = 1.0 x 10~ 7 T[cii] (T 0 ). For the remainder of 


this paper, uncertainties on Lfir (and inferred SFR) only 
include the flux measurement uncertainties, not that of 
the dust temperature and luminosity to SFR conversion. 

The synthesized beam sizes are 0'.'63 by 0'.'45 for 
J0055+0146 and 0"76 by 0"64 for J2229+1457. The bet¬ 
ter resolution for J0055+0146 is mostly due to higher 
elevation of observation. We used the CASA IMFIT task 
to fit 2D gaussian models to these maps. For J0055 both 
















































4 


Willott et al. 


the continuum and line are resolved with deconvolved 
source sizes of O'/51 ±0'/13 by 0'/35±0'/26 at position an¬ 
gle 87 degrees and 0"50±0'/14 by 0'/18±0'/27 at position 
angle 62 degrees, respectively. At the distance to this 
quasar the spatial extent of 0.5" is equal to a linear size 
of 2.9 kpc. We note that the missing data in the red wing 
of the [C ii] line may cause a bias in the size and inclina¬ 
tion if the emission comes from a rotating disk, but there 
is no evidence for this based on the similarity of the line 
and continuum sizes. For J2229+1457 the continuum is 
too poorly detected to attempt a size measurement and 
the line emission is only marginally more extended than 
the beam size. In several other z ~ 6 qu asars velocit y 
gradients across the sources are observed (Wil3 Wang 
et al. 2013). Velocity gradients are not seen for either of 
these two quasars, although for J0055+0146 the missing 
data for 150 km s ^ 1 of the red wing hampers our ability 
to detect such a gradient. 

4. DISCUSSION 

4.1. Evolution of far-IR luminosity 

In Wil3j we reported on the low far-IR luminosities of 
the two previously observed CFHQS quasars and impli¬ 
cations for the relatively low SFR of these quasar host 
galaxies relative to the black hole accretion rate. We 
now revisit this issue with the sample of four z ss 6 
CFHQS quasars with ALMA observations. We note that 
this sample includes four of the six CFHQS quasars with 
measured black hole masses within the absolute mag¬ 
nitude range —25.5 < M 1450 < —24 at a declination 
low enough for ALMA observation. The two unobserved 
quasars have 7 x 10 s < Mbh < 1O 9 M 0 , and were not 
observed due to limited time available and the desire to 
study the lowest mass black holes from CFHQS. There¬ 
fore there is a slight bias to low black hole mass in this 
sample compared to pure UV-luminosity-selection. 

Two of the four quasars are well detected in the contin¬ 
uum with fluxes of 211 ± 34 (J0055+0146) and 120 ± 35 
(J0210-0456) fiJy. J2229 has a marginal 2a detection of 
54±29 /iJy (Figure[3]and Tablefl]) and J2329-0301 is un¬ 
detected with a la rms of 30 /iJy. We combine the four 
values of far-infrared luminosity derived from these mea¬ 
surements assuming that J2329-0301 has a flux equal to 
its 2cr upper limit of 60 /iJy. The mean and standard de¬ 
viation of the sample is Lfir = (2.6± 1.4) x 10 11 A 0 . We 
note this is much lower than the values of 10 12 — 10 13 A 0 
typically discussed for z ss 6 quasars due to two factors, 
firstly that the CFHQS sample here have lower AGN 
luminosity than most known z ~ 6 quasars and a cor- 
relation between AGN and far-IR lumino sities is present 
(|Wang et al.||2011| |Omont et al.||2013|), but also that 
our small sample is selected on quasar rest-frame UV 
luminosity and black hole mass, whereas previous stud¬ 
ies have focussed on quasars with pre-ALMA millimeter 
continuum detections. 

The implication is that these quasars have very high 
black hole accretion rates as inferred from the AGN bolo- 
rnetric luminosity, Lb 0 i, yet relatively low SFR. Such 
a scenario is consistent with the well-known evolution¬ 
ary model whereby the optica l quasar phase comes after 
the main st ar forming phase (Khandai et al.||2012 Lapi 
et al. |2014 ), possibly due to quasar feedback inhibiting 
gas cooling and star formation. The measured ratio of 
Afir/AboI = 0.035 for the four CFHQS quasars is only 



Fig. 4. — Stacked mean far-infrared luminosity for samples of 
quasars at different redshifts. Details of the samples are described 
in the text, but all samples are selected to include roughly the same 
range in bolometric luminosity centred on LboI ~ 7 x 10 12 Lq, the 
mean LboI of the four z « 6 CFHQS quasars plotted with the blue 
square. There is a clear rise in Lpi^ up to a peak at 2 < z < 3 
followed by a decline to z = 6. The magenta c urve shows t he mean 
Ij |r 11 ; due to star formation for the model oflVeale et al. (2014) 
including scaling by a factor of 2 to account for stellar mass loss. 


found in the optical quasar phase of co-evolution at a 
time ~ 1 Gyr after the ons et of activity for the z = 2 
model of Lapi et al. (20141. Given that this is the age 
of the universe at z = lTtlie evolution must occur more 
rapidly at higher redshift. However, the effect of AGN 
variability may also be important leading to a selection 
effect whereby AGN luminosity-selected objects are ob¬ 
served to have l ower ratios of Afir/Abqi than the tim e- 
averaged values ( Hickox et al.| | 2014 |Veale et al.||2014j ). 

We have previously shown (|Wil3p that the two z = 6.4 
CFHQS quasars observed with ALMA in Cycle 0 have 
Afir lower than quasars of similar AGN luminosity at 
lower redshift. On the other hand, at fixed AGN lurni- 
nosity Afir is observed to rise from z = 0 to z = 3 


( Serjeant et al.||2010 Bonfield et al.poTlj Rosario et al. 
20121120131). We next analyze the evolution of Afir us- 


2012| 2013j). We next analyze the evolution of Afir us¬ 
ing our expanded ALMA sample at z = 6 and compa¬ 
rable low-redshift data. At all redshifts we determine 
the mean Afir for optically-selected quasars and X-ray 
AGN in a narrow range of Ab 0 i corresponding to the 
mean Ab 0 i of the four CFHQS z ~ 6 quasars in this 
paper (A Bo i ~7x 10 12 A 0 ). 

At z < 3 we use three datasets b ased on Her- 
schel imaging of AGN. Serjeant et al. (j2010j) stacked 
Herschel SPIRE data of optically-selected quasars and 
quoted their results as rest-frame 100 /.inr luminosity. 
We adopt Afir, = 1.43 vL v (100 /im) ([Chary fe Elbaz| 


20011 to convert to far-infrared luminosity. The abso- 


lute magnitude bin —26 < Tab < —24 corresponds well 
to AboI ~ 7 x 10 12 Lq and we have trimmed the size of 
the highest redshift bin from 2 < z < A to 2 < z < 2.7 



quasars in COSMOS. Due to the shorter wavelength of 
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PACS than SPIRE they presented results in rest-frame 
60 /im luminosity. W e adopt Afir = 1.5 uL u (G0 /tun) 
(Chary & Elbaz 2001) to convert to far-infrared lumi¬ 
nosity From this study we use only the highest red- 
shift, highest lumin osity bin as th is com pares well with 
A Bol ~7x 10 12 Lq. Rosario et al. (2012) determined the 
mean infrared luminosity with PACS for X-ray-selected 
AGN from the COSMOS survey. We note that the X-ray 
selected AGN sample contains a mixture of broad-line, 
narrow-line and lineless AGN and these may have dif- 
ferent evolutionary properties, but Rosario et al. ( [2013 1 
showed that the mean Afir of quasars and X-ray-selected 
are simila r at a given AGN l uminosity and redsliift. We 
use the Rosario et al. (2012) data from the AGN lumi¬ 
nosity bin 8 x 10 n < AroI < 2 x 10 13 A 0 . Whilst most 
sources in this bin have AroI < 3 x 10 12 Aq, we consider 
the results appropriate to compare to the z ~ 6 quasars 
as the correlation bet ween Afir and Ar q i is very shallow 
at this luminosity in Rosario et al. (2012). 

Figure [4] plots data from these three low-redshift stud¬ 
ies with the CFHQS ALMA bin at 6 < z < 6.5. The 
three low-redshift studies show a rise in Afir of a factor 
of 4 from z = 0.3 to z = 2.4. This rise is attributed to the 
general increase in massive gala xy specifi c st ar form ation 
rate over this redshift range (Hickox et al. 2014). The 
6 < z < 6.5 bin has a large dispersion due to the range 
in 1.2 mm continuum flux measured for the 4 quasars. 
The mean Afir at z « 6 is a factor of about 2 lower 
than the z = 0.3 bin and 6 lower than the z = 2.4 peak 
at the so-called quasar epoch. There is clear evidence 
here for a turnaround that mimics the evolution of the 


quasar luminosity function ( McGreer et ah 2013) and 
star formation rate density (Bouwens et al. 2014), albeit 
with a much less steep high-redshift decline due to the 
fact we are measuring star formation in special locations 
within the universe where dark matter halos must have 
collapsed much earlier than typical in order to build up 
the observed black hole masses of ~ 10 s M 0 . 

What is the physical reason for this turnaround at 
z > 3? In the evolutionary picture where the opti¬ 
cal quasar phase follows the starburst phase one would 
expect the star formation and black hole accretion to 
be more tightly coupled at high-redshift where there is 
barely enough time for star formation to have decreased 
substantially. A clue may come from one of the few differ- 
ences b etween quasars at these two epochs. Willott et al.| 
|2010a) showed that the Eddington ratios of matched 
quasar luminosity samples at z = 2 and z = 6 are signifi¬ 
cantly different with the z = 6 quasars having a factor of 
3x higher Eddington ratios and therefore 3x lower black 
hole masses than at z = 2. Such a difference exists be¬ 
tween the typical black hole mass of our CFHQS ALMA 


sample and that of the highest luminosity bin of Rosario 


et al. (2013). T his Eddington ratio evolution is observed 
in other studies (De Rosa et al.|20TT Trakhtenbrot et al. 


2011 |Shen fe K elly 201 (2) and predicted by many theo¬ 
retical works due to t he increase in gas s upply to black 
holes at high-redshift (|Sijacki et al.]|2014|). 

In Figure [4] we also plot a theoretical curve of mean 
Afir versus redshift for a simulated sample of rest-frame 
UV-selected quasars in the same Ar 0 i range as the ob - 
served quasar samples for the model of Veale et al.[ 
(2014). This model assumes an evolving linear relation- 
ship between star formation and black hole growth. The 


variant of the model plotted here is the “accretion” model 
where the quasar luminosity is proportional to the black 
hole growth rate and the Eddington ratio distribution is a 
truncated p ower-law with slope /3 = 0.6 (dashed curve in 
Figure 8 of Veale et al. 2014). The model is constrained 
by the observed evolving quasar luminosity function and 
the local ratio of black hole to galaxy mass. We have 
scaled this model with a factor of 2 x increase in Afir to 
account for stellar mass loss. 

As seen in Figure [4] this curve increases from low red¬ 
shift to the peak quasar epoch at 2 < 2 < 3 by about the 
same factor as the data, although the total normalization 
of the curve is lower by a factor of 3 to 4. [Veal e et al.| 
(20141 discuss some of the reasons why the normaliza¬ 


tion may be lower than the observations. The decrease 
in Afir with increasing cosmic time from z = 2 to z = 0 
for fixed luminosity quasars is due to the fact that such 
quasars are rarer at lower redshift and on the steep end 
of the luminosity function where scatter is more impor¬ 
tant. This behaviour also follows from the general de¬ 
crease in specific star formation rate with cosmic time. 
The high-redshift behaviour of a decline from z = 3 to 
z = 6 matches our observations, so it is instructive to 
understand why this occurs in the model. It is due to 
the assumed (1 + z) 2 evolution of the ratio of accretion 
growth to stellar mass growth, but this assumed evolu¬ 
tion is also degenerate with evolution in the Eddington 
ratio. As discussed previously there is observational evi¬ 
dence for positive evolution in the Eddington ratio from 
z = 2 to z = 6, meaning that the ratio of accretion 
growth to stellar mass growth may change more gradu¬ 
ally than (1 + z) 2 . 

A possible alternative explanation for the low Afir at 
z = 6 is that at these early epochs insufficient dust has 
been generated so that star form ation occurs more often 


et al. 2013 

Fisher et al. 

2014 

Uta et af. 2014) 


In this 


SFR of a luminous quasar hosting galaxy. However, two 
lines of evidence point towards this not being the main 
factor for our quasar sample. First, quasars at z = 6 
are known to have emission line ratios similar to lower 
redshift quasars inferring hi gh metallicity at least close 


to the accreting black hole (jFreudling et al. 2003). Sec¬ 
ond, the [Cn] luminosities in three of the four quasars 
are high (see below), suggesting high carbon abundances 
throughout the host galaxies. 

4.2. The [C II/ - far-IR luminosity relation 

Three of the four CFHQS ALMA quasars are detected 
in both [C Ii] line and 1.2 mm continuum emission (tw o 
new detections in this paper and J0210—0456 in Wil3). 


With the low Afir discussed in the previous section, 
these quasar host galaxies probe a new regime in Afir 
at high-redshift. In Figure [5] we plot the ratio of [Cn] 
to far-IR luminosity as a function of Afir . Also plot¬ 
ted are several samples from the literature which, due 
to ALMA at high-redshift and Herschel at low-redshift, 
are rapidly increasing in size and data quality. The low- 
redshift z < 0.4 sample of galaxies is from |Gracia-Carpio| 
et al. (2011 and in prep.) and contains a mix of normal 
galaxies, starbursts and ultra-luminous infrared galaxies 
(ULIRGs), some of which contain AGN. The ULIRGs 
show a [Cll] deficit that has been widely discussed in 
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Fig. 5.— The ratio of [Cn] to far-IR luminosity as a function 
of far-IR luminosity. High-redshift (z > 5) sources, mostly quasar 
host galaxies, are identified with large symbols. Error bars are only 
plotted for the CFHQS ALMA sources to enhance the clarity of the 
figure. The solid and dotted blue lines show the best fit power-law 
and la uncertainty for the z > 5 sources. The z > 5 relation is 
largely consistent with the distribution of data at lower redshift. 


contamination of Lftr (Sargsyan et al. 

2012 ), high gas 

fractions 

Gracia-Carpio et al. 

2011 ) or the dustiness ? 

temperature and/or density of star forming regions (Far- 

rah et al. 

2013 Magdis et al. 2014). P 

revious observa- 


1 , 1011 s oi mgn i^FTR z > 5 SUSS quasars (Maiolino et al. 
2005 Wang et al.||2013) showed a similar deficit. How- 


ever many 0.5 < z < 5 ULIRGs do not show this deficit 
and have L[cii]/Tfir ratios comparable to low-redshift 
star-forming galaxies (Stacey et al.|20f0). Thi s is visible 


in Figure [5] for the 0.5 < z < 5 compilation of |De Looze| 


et al. (20141). 

By adding three z > 6 quasars with 10 11 < Tfir < 
10 12 l q to Figure [5] we have greatly expanded the range 
of luminosities at the highest redshift . Large symbols on 


Looze et al. 

2014D sources are 

HLSJ091828.6+514223 

at 2 = 5.24 

Rawfe et al. 

2014 

) , HFLS3 at 2 = 6.34 

(Riechers et a 

. 2013) and 

DSS J1148+5251 at 2 = 6.42 


mostly selected tor followup based on high L F ir. The 
quasar ULAS J1120+0641 at z = 7.1 has a more moder¬ 
ate I/fir and -Qciii/Tfir ratio and lies in between the 
CFHQS and high Lfir objects on the plot. 

Although we are wary of the selection effects in Fig¬ 
ure [5] and the as yet unknown cause for the change in 
L[cii]/Lfir with Lfir at high-redshift, the new data pro¬ 
vide the opportunity to make the first measurement of 
the slope of this relation at z > 5. We fit the 12 z > 5 
sources with a single power-law model of the dependence 
of Tfir on L[cn] incorporating the observational (but 
not systematic) uncertainties using a MCMC procedure. 
The best-fit relation is 

logio ^fir = 0.59 + 1.271og 10 L[cii]- (1) 

Therefore the L[cii]/Tfir ratio line plotted in Figure [5] 
has a logarithmic slope of (1/1.27) — 1 = —0.21. The 
MCMC la uncertainties, based solely on the observa¬ 


tional data, are plotted as dotted lines. These also favor 
a shallow negative slope, not nearly as steep as the slope 
that would be fit to the previous 2 > 0.5 data that cov¬ 
ers only a narrow range of Tfir- The L[cii]/Afir ratios 
of the CFHQS and ULAS 2 > 6 quasars are not greatly 
d ifferent to those of similar Lfir galaxies at low-redshift. 


Wang et al. (2013) note that the low L[ Cii i/Lfir ratios 
for SUSS 2 > 0 quasars may be at least in part due to 
AGN contamination of the far-IR emission. Future ob¬ 
servations at higher spatial resolution will be critical to 
examine differences in the spatial distr ibution of the line 
and continuum emission (e.g. Cicone et al. 2014). We 


expect to observe that the dust continuum is more com¬ 
pact than the [C 11 ] line in the high Lfir 2 > 5 quasars, 
due to either more centrally-concentrated starbursts with 
higher dust temperatures, like local ULIRGs, or AGN 
dust-heating. In contrast we expect the low Lfir 2 > 5 
quasars have star formation spread more evenly through¬ 
out their host galaxies, with similar spatial distribution 
of line and continuum emission. 

4.3. The 2 ~ 6 .Mbh ~ cr and Mbh — Md yn relationships 

The combination of black hole mass estimates and [C 11 ] 
line host galaxy dynamics for these 2 « 6 quasars allows 
us to investigate the black hole - galaxy mass correlation 
at an early epoch in the universe. The evolution of this 
relationship is a critical constraint on the co-evolution 
(or not) of galaxies and their nuclear black holes. As 
discussed in the Introduction there are reasons to be¬ 
lieve that past studies using only t he most m assive black 


holes from SDSS quasars (e.g. Wang et al. 2010) were 
prone to a bias where one would expect the black holes 


et al. 

2005 

Lauer et al. 

2007) 


as observed. With new 


hypothesis and determine any real offset from the local 
relationship. Additionally, most previous work in this 
area has used the molecular CO line to trace the gas dy¬ 
namics. CO is usually more centrally concentrated than 
[C 11 ], so [C II] potentially probes a l arger f racti on of the 


total mass (although we note that Wang et al. 2013 
found similar dynamical masses using CO and [0 11 ] for 
their 2 « 6 quasars). 


In addition to the CFHQS and Wang et al. (20131 
quasars we add to our study two other 2 > 6 quasars 
observed in the [Cli] line: SDSS J1148+5251 at 2 = 6.42 
and the most distant known quasar, ULAS J1120+0641 
at 2 = 7.08. Both these quasars have Mg 11 -derived black 


hole masses (De Rosa et al. 2011 2014) with very low 


measurement uncertainties. The black hole masses for 
all thre e CFHQS quasars als o come from Mg 11 measure¬ 


ments (Willott et al. 1 2010a). Some of these spectra are 
of moderate SJNR ana have substantial measurement un¬ 
certainties on the black hole masses. To all the quasars 
with Mgn-derived black hole masses we add a 0.3 dex 
uncertainty to the measurement uncertainties to account 
for the dis pers ion in t he rev erberation-mapped quasar 
calibration ( Shenet_aL|2008 k_ 


None of the Wang et al. 7 2013) quasars have Mg II mea¬ 


surements so black hole masses are estimated assuming 


served for most 2 ~ 6 quasars 

((Jiang et al. 2007 

Kurk 

et al. 2007 Willott et al. 2010a 

IJ^Rosa eHil.] 2011 

The 
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at z « 6 is 0.3 dex (Will ott et al.[201 0a). We add 0.3 dex 
uncertainty from the observed dispersion in the Edding¬ 
ton ratio distribution in quadrature to the 0.3 dex due to 
the dispersion in the reverberation-m apped quasa r cali- 
bration for a total uncertainty on the Wang et al. (20131 
quasar black hole masses of 0.45 dex. 

First we consider the Mbh - o relationship. For nearby 
galaxies a is the velocity dispersion of th e galaxy bulge. 


At high-redshift bulges are less common (Cassata et al. 
2011) and we do not expect the z ~ 6 kinematics to 


match that of a pressure supported bulge. With the lim¬ 
ited spatial resolution of current data we cannot be sure 
the [C ii] gas is distributed in a rotating di sk, alt hough 


there is evi d ence o f this for some sources (Wang et al. 
2013). Ho] (2007a) discusses the relationship and cali¬ 
bration of bulge velocity dispersion and disk circular ve¬ 
locity and concludes that although there is additional 
scatter one can relate molecular or atomic gas in a disk 
to stellar bulges. The major complication is the incli¬ 
nation of the disk. For a random sample of inclinations 
this can be modelled, however there is a possibility that 
quasars have disks oriented more often fa ce-on , reducing 
the line-of-sight velocity dispersion (|Ho| 2007b). 

We determine a using the metlio< 


(2007b), 


specifically setting the [C n] line full-width at 20% equal 
to 1.5 x the FWHM as expected for a Gaussian since most 
of the lines are approximately Gaussian. The [C Ii] emit¬ 
ting gas is assumed to be in an inclined disk where the 
inclination angle, i. is determined by the ratio of minor 
(flmin) and major (a maj ) axes, i = cos _1 (a min /a maj ). The 
circular velocity is therefore u c j r = 0.75FWHM[ CII ]/sin*. 
For all of the quasars in this study we determine an in¬ 
clination from the [C n] data or assume an inclination if 
one or both th e major an d minor axes are unresolved. All 
the quasars in Wang et al. (2013) were spatially resolved, 
although some had quite large uncertainties on a m i n and 
a ma j. We adopt the inclination angles from their paper. 

For the [C Ii] emission of J0210—0456, a ma j = 0752 ± 
0"25 (2.9 kpc) with i = 64° (fWil3). For J0055+0146, 


the quasars are distributed around the local relationship 
rather than all being offset to low a as is co mmonly be¬ 
lieved to be the cas e. As noted by Wang et al. (2010), 
using the method of Ho (2007b), rather than calculating 
a as FWHM/2.35, leads to much higher a. Note that 
this is without adopting extreme face-on inclinations for 
most quasars. There are still several quasars, such as 
SDSS J1148+5251 and ULAS J1120+0641, that have val¬ 
ues of a considerably lower than the local relation. 

The main result of Figure [6] is that whilst there is lit¬ 
tle mean shift between the z = 0 and 2 « 6 data, there 
is a much larger scatter in the data at z ~ 6, well be¬ 
yond the size of the error bars. This larger scatter at an 
early epoch is expected based on dynamical evolution, 
incoherence in AGN/starburst activity and t he tighten - 
ing of the relation over time from merging (Peng [2007). 
We note that our hypothesis that the bias described in 
the Introduction would lead to the lower Mbh quasars 
being located on the local relation with a lower scatter 
than the high Mbh quasars is not supported by these 
observations. The scatter in log 10 a at Mbh ~ 10 s M 0 
is about the same as that at Mbh > 10 9 M 0 

We go one step further from a to determine dynam¬ 
ical masses using the deconvolved [Cii ] sizes. For con - 
sistency, we follow the method of Wang et al. (2013). 
The dynamical mass within the disk radius is given by 


M 


dyn 


1.16 x 10 5 u 9 ir D Mq where D is the disk diame¬ 
ter in kpc and calculated as 1.5x the deconvolved Gaus¬ 
sian spatial FWHM. The resulting dynamical masses are 
given in Table [2] We note that there is considerable 
uncertainty on these values due to the unknown spatial 
and velocity structure of the gas, the marginal spatial 
resolution and limited sensitivity that means we may be 
missing more extended gas. Due to the these uncertain¬ 
ties we do not pl ace formal err or bars on the dynamical 
masses, following Wang et al. (2013). Higher resolution 
data in the future are required to confirm the derived 
masses. 

SDSS J1148+5251 has been extensively studied in [C nj 


Umaj = 0750 ± 0714 (2.9 kpc)i with i = 69° (Section 3, (Maiohno_j£jLL|j2005; Walter et al. 2009 Maiolino et al. 


assuming no bias from missing red wing data). The 
[C ii] emission of J2229+1457 is only marginally spa¬ 
tially resolved (0785 versus beam size of 0776) and we 
estimate an intrinsic FWHM of « 074 (2.4kpc). An in¬ 
clination of i = 55° is assumed as this is the median 
incli natio n angle for the resolved sources in this paper 
and|Wang et al.| (|2013|). Neither SDSS J1148+5251 nor 
I I.AS.) 1120 • OO'll have published inclination angles, so 
we also assume i = 55° for both of them. We adopt 
FWHMfcm - 287 ± 28kins" 1 for SDSS J1148+5251 
(Walter et al. 2009) an d FWHM[ CII | = 2 35 ± 35kms _1 
tor ULAS J1120+0641 (Venemans et al. 2012) . Values 
of black hole masses and a tor this sample are provided 
in Tabled 

Figure [6] shows the Mbh — cr relationship for the z ~ 6 
quasar sample. Uncertainties on black hole masses in¬ 
clude the scatter in the calibration as described previ¬ 
ously. Uncertainties in a include FWHM measurement 
uncertainty plus a 10% uncertainty for the conversion 
from FWHM to u c ; r and 15 % for the co nversion from *) c j r 
to a as seen in the sam ple of Ho (2007a). The black line is 
the local correlation of Kormendy & Ho (2013) with the 


2012| Gicone _et al. 2014+ The highest resolution ob¬ 
servations by Walter et al. (2009) revealed a very com¬ 
pact circumnuciear star burst with radius 0.75 kpc and 
FWHM[cn] = 287 ± 28kms~ 1 . For an assumed inclina¬ 
tion of i = 55° this give s _Mdy n = 1.8 x 10 10 M 0 . For 


comparison Walter et al. (2004) determined a dynamical 
mass from CO emission in this quasar of 5.5 x 10 10 M 0 
within a larger radius of 2.5 kpc, the larger radius being 
the main difference between the results. Recent observa¬ 
tions have shown more complex [C n] emission including 
evidence for gas extended over tens of kpc and at high 


veloci ties in dicative of outflow (Maiolino et al. 


cone et al. 


2012 Ci- 

2014b. We adopt M dyn = 1.8 x lO 10 M, 


lor 

SDSS J1148+5251 noting that the true value could be 
several times larger. The most distant known quasar, 
ULAS J1120+0641 at 2 = 7.08, has been well detected 


in [Cii] , alt hough not yet spatially resolved (Venemans 
et al. 2012). Based on the published FWHM[ CII ] = 


gray band the ±lcr scatter. The first thing to note is that 


235 ± 35kms _1 and assuming a spatial FWHM of 3 kpc 
(similar to the other quasars resolved by ALMA) and 
i = 55° we determine M dyn = 2.4 x 10 10 M 0 . 

In Figure [7] we plot black hole mass versus galaxy dy¬ 
namical mass for the most distant known quasars. The 
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TABLE 2 

Mass determinations for z > 5.7 quasars 


Name 

2 [CII] 

M bh (M 0 ) a 

a b 

Mdyn (Mg) 

Refs. c 

CFHQS J0055+0146 

6.0060 ± 0.0008 

(2-4 HI) 

X 

o 

00 

207 ± 45 

4.2 x 10 10 

1,2 

CFHQS J0210—0456 

6.4323 ± 0.0005 

(0-8±S;g) 

00 

o 

X 

98 ± 20 

1.3 x 10 10 

1,2,3 

CFHQS J2229+1457 

6.1517 ± 0.0005 

(i-2tS; 4 8 ) 

X 

o 

00 

241 ± 51 

4.4 x 10 10 

1,2 

SDSS J0129-0035 

5.7787 ± 0.0001 


x 10 s 

112 ±21 

1.3 x 10 10 

4 

SDSS J1044-0125 

5.7847 ± 0.0007 

(1-lto.?) 

O 

o 

X 

291 ± 76 

__ d 

4 

SDSS J1148+5251 

6.4189 ±0.0006 

(4-9±|;§) 

x 10 9 

186 ± 38 

1.8 x 10 10 

5,6,7 

SDSS J2054-0005 

6.0391 ± 0.0001 

(0.9±S:S) 

x 10 9 

364 ± 67 

7.2 x 10 10 

4 

SDSS J2310+1855 

6.0031 ± 0.0002 

(2.81®*) 

x 10 9 

325 ± 61 

9.6 x 10 10 

4 

ULAS J1120+0641 

7.0842 ± 0.0004 

(2-4±?;|) 

x 10 9 

144 ± 34 

2.4 x 10 10 

8,9 

ULAS J1319+0950 

6.1330 ±0.0007 

(2.1±f|) 

X 

o 

CO 

381 ± 91 

12.5 x 10 10 

4 


Notes.— 

a Derived from Mgn A2799 observations if possible, else from Eddington luminosity assumption. Uncertainties include observational errors 
plus systematics based on calibrations. 

b Derived from Gaussian FWHM fit to [C ii] spectrum using method of|Ho|(|2007b|) including an inclination correction (see text for individual 
inclinations assumed). Uncertainties include observational errors plus systematics based on calibrations. 

(6) 


References: (1) T his paper, ( 2) Willo tt et a l. 
| Walter et al.| ( |2009) , (7) |De Rosa et~aT] ( |2Ulf| 


I] (|2010ab 

Tyo) | Vene: 


(3) iWillott et al . 

Venemans et al. |2(^^ ~ 


ld20 13|), (4) | Wang et alT fl2013|), (5) |Maiolino et al. | (|2005|), 
(9) |De Rosa et ak|fl2UU ). 


This quasar does not have a dynamical mass calculation in Wang et al. (2013) due to the difference in the [C n] and CO line profiles. 




Fig. 6. — Black hole mass versus veloci ty d i spersion calculated 
from the [C n] line using the method of |Ho| ([2007b) for z « 6 
quasars. Quasars from the CFHQS are shown as blue squares and 
the other symbols show quasars from the SDSS and ULAS surveys. 
The black line with gray shading is the local correl ation zb la scatter 
of bla ck hole mass and bulge velocity dispersion ([Kormendy &; Ho| 
|2013[ ). The z « 6 quasars are distributed around the the local 
relationship, but with a much larger scatter and some quasars with 
significantly lower a for their Mbh- 


black line and gray shading represent the local correla- 


tion o f Mbh with bulge mass Mb u i ge (Kormendy & Ho 


2013). In the absence of gas accretion and mergers the 


present stellar bulge mass represents the sum of the gas 
and stellar mass at high-redshift, so it is a good com- 
paris on for the dynamical m ass within the central few 
kpc. Kormendy & Ho (2013) note that their correlation 
(their equation 10] gives a black hole to bulge mass ra¬ 
tio of 0.5% at Mbuige = 10 11 M 0 that is 2 to 4 times 
higher than previous estimates due to the omission of 


Fig. 7. — Black hole mass versus host galaxy dynamical mass 
for z £3 6 quasars. Symbols as for Figure [m The black line with 
gray shading is th e loca l correlation drier scatter from the work of 
|Kormendy & Ho| (|2013) equating M^ yn to Mb u i ge - The CFHQS 
quasars he on the local relationship and do not show the large offset 
displayed by the most massive black holes. Uncertainties in Md yn 
have not been calculated due to the reasons given in the text. 


pseudobulges, galaxies with uncertain A/bh and ongoing 
mergers. 

The position of the high -2 data with respect to low 
redshift is fairly similar to Figure [6j not surprising be¬ 
cause v C i r derived from the [Cll] velocity FWHM is a 
major factor in both er and M ( j yn . The points are shifted 
somewhat further from the local bulge mass than for the 
local velocity dispersion. This shift is due to the smaller 
size of galaxies at high-redshift, as the size is the only 
term in the derivation of dynamical mass not in a. We 
note the much greater dynamic range in black hole mass 
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(2 dex) than in dynamical mass (1 dex) in our sample. 
This is likely due more to our selection over a wide range 
of quasar luminosity than to a non-linear relationship 
between these quantities at z = 6. 

All three of the CFHQS quasars lie within the l ocal l cr 
scatte r and the one M BH ~ 10 8 M 0 quasar in Wang et al. 
(2013) is only a factor of 4 greater than the local rela¬ 
tionship. In contrast the Mbh > 10 9 M 0 quasars tend to 
show a larger scatter and large r offset above the local re¬ 
lationship as previously found (Walter et al. 2004 Wang 
et al. 2010; Venemans et al. 2012; Wang et al. |2013|). We 


caution that there are considerable uncertainties m some 

of these measurements as already discussed, but in dy¬ 
namical mass the results look more like we would expect 
based on the quasar selection bias effect. 


5. CONCLUSIONS 

During ALMA Early Science cycles 0 and 1 we have 
observed a complete sample of four z > 6 moderate lumi¬ 
nosity CFHQS quasars with black hole masses ~ 10 s M 0 . 
Three of the four are detected in both far-IR continuum 
and the [Cn] emission line. The far-IR luminosity is 
found to be substantially lower than that of similar lumi¬ 
nosity quasars at 1 < z < 3. Assuming that far-IR lumi¬ 
nosity traces star formation equally effectively at these 
redshifts this implies that at z ~ 6 quasars are grow¬ 
ing their black holes more rapidly than their stellar mass 
compared to at the peak of the quasar epoch (1 < z < 3). 

The ratios of [CII] to far-IR luminosities for the 
CFHQS quasars lie in the range 0.001 to 0.01, similar 
to that of low-redshift galaxies at the same far-IR lumi¬ 
nosity. This suggests a similar mode of star-formation 
spread throughout the host galaxy (rather than in dense 
circumnuclear starburst regions that have lower values 
for this ratio in local ULIRGs). Combining with pre¬ 
vious 2 > 5.7 quasar data at higher Lfir we find that 
the far-IR luminosity dependence of the [Cll]/FIR ra¬ 
tio has a shallow negative slope, possibly due in part to 
an increase in Lfir due to quasar-heated dust in some 
optically-luminous high -2 quasars. 

The three CFHQS quasars well-detected in the [Cn] 
emission line allow this atomic gas to be used as a tracer 
of the host galaxy dynamics. Combining with published 
data on higher black hole mass quasars we have investi¬ 
gated the Mbh — cr and Mbh — Md yn relations at 2 w 6. 
We show that the 2 = 6 quasars display a Mbh — u rela¬ 


tion with similar slope and normalization to locally, but 
with much greater scatter. Similar results are obtained 
for the Mbh - M dy n relation with a somewhat higher nor¬ 
malization at 2 = 6 and a higher scatter at high Mbh- 
As discussed in 

Combining our results on the relatively low Tfir. for 
Mbh ~ 10 8 M 0 z ss 6 quasars with their location on 
the Mbh — cr relation leads to something of a paradox. 
The fact these quasars lie on the local Mbh — a relation 
suggests that their host galaxies have undergone consid¬ 
erable evolution to acquire such a high dynamical mass. 
So why is it that this mass accumulation is not leading 


to a high star formation rat e? As discussed in Wil3 
simulations such as those of Khandai et al. (2012|) and 
Lapi et al. (2014) predict that such low ratios of SFR to 
black hole accretion occur after episodes of strong feed¬ 
back that inhibits star formation throughout quasar host 
galaxies. Another possibility mentioned in Section 4.1 is 
that Lfir fails to trace star formation so effectively in 
these hi gh-r e dshif t galaxies, due to lower dust content 
(e.g. Ouchi et al. 2013). Note that using T[cn] as a star 
formation rate tracer instead of Lfir, would give higher 
SFR by a factor of three for one of the CFHQS quasars. 

Higher resolution follow-up [C ii] observations of these 
quasars are critical to measure more accurately the dis¬ 
tribution and kinematics of the gas used as a dynamical 
tracer in order to reliably determine the location and 
scatter of the correlations between black holes and their 
host galaxies at high-redshift. 


Thanks to staff at the North America ALMA 
Regional Center for processing the ALMA data. 
Thanks to Melanie Veale for useful discussion and 
providing her models in electronic form. This 
paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2012.1.00676.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF 
(USA) and NINS (Japan), together with NRC (Canada) 
and NSC and ASIAA (Taiwan), in cooperation with the 
Republic of Chile. The Joint ALMA Observatory is oper¬ 
ated by ESO, AUI/NRAO and NAOJ. The National Ra¬ 
dio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. 

Facility: ALMA. 
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